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Dear Sir: 

We; hereby declare that: 

(1) . We, 

James V. Candy, 417 Brigham Lane, Danville, CA 94526, a citizen of the 

United States; 

and 

David H. Chambers, 1665 Altamar Way, Livermore, CA 94551, a citizen of 
the United States; 

are citizens of the countries and residents of the cities and states 
identified. 

(2) . We are the inventors named in the subject application. 
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Declaration of Prior Invention to Oveixoine Cited Reference 



(3) . The subject application SN 10/661,249 was filed September 11, 2003. 
Tlie subject application claims benefit of provisional patent application SN 
60/410,575 filed September 12, 2002. 

(4) . In the Office Action mailed October 2, 2007, the claims in the subject 
application were rejected over references that include the Kerbrat et al reference 
(Transactions of Ultrasonics, August 2002). The date of the Kerbrat et al reference 
is on or about August 1, 2002. Accordingly, August 1, 2002 is the date that must 
be overcome. 

Overcomin g the August 1, 2002 Date of Kerbrat et al Reference 

(5) . We made the invention described and claimed in the subject patent 
application (The Invention) in this country prior to the August 1, 2002 date of the 
Kerbrat et al Reference. We made written descriptions of The Invention, we 
made drawings of The Invention, we made reports about The Invention, we 
made tests of The Invention, and we disclosed The Invention to others; all of the 
foregoing were done in this country prior to August 1, 2002. 

(6) . We conceived The Invention in this country prior to August 1, 2002. 
Attached are documents supporting our statements that we conceived The 
Invention in this country prior to the August 1, 2002 date of the Kerbrat et al 
reference. 

RECORD OF INVENTION signed October 18, 2001 

(7) . The attached October 18, 2001 RECORD OF INVENTION was 
completed and signed by us October 18, 2001 well before the August 1, 2002 date 
of the Kerbrat et al Reference. Our signatures, witness's signatures, and the dates 
of signing appear on page 5 of the RECORD OF INVENTION. 

(8) . The attached October 18, 2001 RECORD OF INVENTION was sent to 
the Lawrence Livermore national Laboratory Office of Laboratory Counsel 
(Intellectual Property Law Group I.P.L.G.) and received there October 22, 2001. 
The "RECEIVED" stamp with the date OCT 22, 2001 at the top of page 1 shows 
that the RECORD OF INVENTION was received October 22, 2001. 

(9) . The attached October 18, 2001 RECORD OF INVENTION in the 
Conception of the Invention Section XI on page 4, contains an entry for the 
"Conception Date." This date of April 1, 1996 is prior to August 1, 2002. 
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Declaration of Prior Invention to Overcome Cited Reference 



(10) . The attached October 18, 2001 RECORD OF INVENTION in the 
Conception of the Invention Section XI on page 4, contains an entry for the 
"Conception Place/' and the entry of "Lawrence Livermore National Laboratory" 
is in the United States. 

(11) . The attached October 18, 2001 RECORD OF INVENTION in the 
Conception of the Invention Section XI on page 4, contains an entry for the "First 
Sketch or Drawing Date." This date of 7/15/01 is prior to August 1, 2002. 

(12) . The attached October 18, 2001 RECORD OF INVENTION in the 
Conception of the Invention Section XI on page 4, contains an entry for the "First 

Written Description Date." This date of 7/15/01 is prior to August 1, 2002. 

(13) . The attached October 18, 2001 RECORD OF INVENTION in Section 
XV on page 5 contains blanks wherein the Name and Title and Signature of an 
Authorized Derivative Classifier appear and the signature date of October 26, 
2001 is prior to August 1, 2002. 

(14) . We were reasonably diligent in reducing the invention to practice 
throughout the 42 day period between the August 1, 2002 date of the Kerbrat et 
al reference and September 12, 2002 when our provisional patent application SN 
60/410,575 was filed (hereinafter "The Time Period"). 

Activity During the Time Period 

(15) . We continuously worked on testing, developing, and patenting The 
Invention during The Time Period from the Kerbrat et al Reference date of 
August 1, 2002 until September 12, 2002 when our Provisional Patent Application 
was filed. 

Supporting Documents 

(16) . Attached are documents that support our statements that we 

continuously worked on testing, developing, and patenting The Invention 
during The Time Period from the Kerbrat et al Reference date of August 1, 2002 
until September 12, 2002 when our Provisional Patent Application was filed. We 
were continuing to do research and development relating to The Invention 
during The Time Period from August 1, 2002 until September 12, 2002. During 
The Time Period the inventor David Chambers applied for a "Multi-location 
Appointment Extension for LLNL Employee" to work with the inventor James 
Candy who was already working both at LLNL and UCSB. The attached memo 
from Associate Director to Human Resources (Subject: Multi-location 
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Appo nfnu nt f xton'-ion foi LLNL Employee) statc^ " I hi'^ menic^rtindum 
iequeti*-s appiox al foi Dr HcU \d ) i Chambers to extend his assignment a"? 
Resoarth Frif^uioei at I e Dc ^ itn t„nt of L lectronitb ana Comp'Uor 1 uirtu n:^ 
(FQ ), LnAoisitx' oi Cah-o n o bouta -^^^ib^iia (tC^B; loi be Iutu> t ^ i o 
8'l/02- 11/15/02 Daxid Chimhus wil! torfmue to pcUosm thi, . ^Mgmmnt at 
25% time at a salary ot 8365 per moiith. . . , Dr. Chambers is currently 
pel touumg computer s uiatuj-i'-o ^ (.cmbtK p»optU,t i oi jn human I imm i 
suppoit ot a UCSB proiect Pu < >; f .niation of this \s ork is anticipated to take 
2^% hme until the end ot t! e pi c HKi I his wot k ^\ ill reduce time available 
for LLNL projects to 75%." 

XU\ o ^ssotla^oP t\ oi t ^ IK'man Ris;.t iv. s 

M.emQ Aug 1, 20U2 Wri$i;ht to Nelson 21). 

Memo Aug 1, 20U2 Wright to Not chemistry Chambers. 

Memo \u^^, 2()02 R kps \onc io \ot chemisti\ Chambers 

Memo Au.p, 5, 2002 Wright to cron. 

Memo Aug 5, 2002 Wright to Nelson 20. 

Memo ^ept i 20t 2 K ni oiid to C ha. nix i ^ 

Memo Sept 12, 2002 Nelson to Wnght. 



{17). We do not know and do not believe that the invention has been in 
public, ufie or on sale in this country,, or patented or described in a printed 
pu.blicaticvn in this or any i'oreigr; ccusntry tor more than one year prior to our 
application, and we have never abandoned our invention. 

(18). We further declare that all statements made herein of our own 
knowledge arc true and that all statements made on information and belief are 
believed to be true; and further that these statements were made with the 
knowledge that willful false stalements and the like so made are punishable by 
fine or imp,rii>oiuru>T^t, (.>r both, under Section 1001 of Title 18 of the United States 
Code, and. that such willful false statements may jeopardize the validity of the 
application or any patent issued thereon. 

. .:?^RC-,. ^' ^.^t'Cs^t-.-..^ 

lanuary'^;.^008 •■ (Signatiij^e} Declarant: James \C€!andy 



...j.;.^. — 

January. :^._/2008 - (Signature) Declarant: Davui IL Chambers 
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HI. Abstract , 

The ability to noninvasively focus acoustical energy in tissue and directly on tissue masses (tumors, cysts, etc.) 
is the primary function of this invention. The objective is to provide the capabiiity of focusing acoustic energy 
at a de^red location for the purpose of treating tissue mass while minimizing the cdlatwal damage In Sie 
surrounding tissue. This invention will open new frontiers vwth the implication of noninvasive treatment of 
masses in the medical area along with the expanding tedinology of acoustic surgery. 



IV. Uses of the Invention 

List past uses, current uses and potential uses for your invention 

LLNL or Government uses or possibilities for use: 

Tissue mass removal, non-invasive tumor/cyst destruction, acoustic surgery, mass imaging, nondestructive 
evaluation of materials, secure communications, seismic detection of underground structures 
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Background 

Background of the invention, including teclinical problems addressed by it: 

With the advent of high-speed digitizers, ultrafast computers, inexpensive memory, and the ability to 
construct dense acoustic arrays, the feasibility of noninvasive techniques of acoustic surgery offers a tantaiizing 
alternative to current invasive techniques. In its simplest form, the focusing of acoustic energy to destructively 
treat a mass in surrounding tissue appears to be a reasonable approach to noninvasive surgery especially if the 
medium is homogeneous and therefore can be characterized by attenuation and time delays. Thus, in this case, 
it is a simple matter to focus energy at a desired point in space. When the medium is inhomogeneous focusing 
at a desired focal point is more difficult unless some knowledge of the medium exists a-priori. There are two 
basic approaches that we discuss in this invention to focus energy in an inhomogeneous medium: (1) time 
reversal techniques, and (2) triodel-based techniques, and their combination. The time reversal approach 
requires no model in contrast to the model-based approach. We discuss both in this disclosure. 

This Invention discusses the capability of focusing acoustic energy at a desired location for the purpose 
of treating tissue mass while minimizing the collateral damage in the surrounding tissue. Our approach is 
summarized simply in Fig. 1. First, we must first detect the presence of a tissue mass applying acoustic energy 
propagated into the tissue using an array of ultrasonic transducers. The amount of energy scattered by the mass 
depends on its acoustic parameters (density, sound speed, attenuation, etc.). Once it is detected, the mass must 
be localized to detennine its position within the tissue medium. There are a vanety of methods that can be used 
to perform this operation, but each make some underlying assumption about the characteristics of the medium 
(homogeneous, inhomogeneities with a homogeneous medium, etc.) leading to uncertainties, Once detected 

localized, temporal signatures must be developed to "drive" the array and focus increased energy back onto 
the mass through the medium. The increased energy generates heat, which essentially "cooks" the mass 
insuring its destruction. 

We demonstrate that confined focusing in tissue is possible and discuss various approaches to achieve 
this focusing. Through "global" focusing, that is, insonifying a iarge region in the tissue the mass is detected and 
localized, "zonal" focusing is perfomned to extract or "zoom in" on the tissue mass under scrutiny. After it is 
decided to treat the mass, increased acoustic energy is transmitted back onto the mass. 




Figure 1, Conc^tual Ultrasonic Focu^ng System for Noninva^ve Mass Treatment: (1) Scatteiing 
tissue medium (3 large masses shown). (2) Focusing array. (3) Mass treatment via ultrasonic 
focusing. 



Time Reversal Focusing 



When a source propagates through a spatio-temporal medium, the resulting wave front is distorted. If 
the medium is homogeneous and the source resides in the near field, then a spherical-type wave front evolves. 
But if the medium is inhomc^eneous, then a distorted wave front results. In the first case, simple time-delay 
processing is sufficient to enhance the field at a given point; however, for inhomogeneous media the required 
time delays and amplitude are more difficuit to estimate. The use of delay estimation and even adaptive delay 
estimation techniques become quite limited and unsuccessful in an inhomogeneous medium excited by a 
broadband incident field requiring an alternative approach to solve the focusing problem. A viable alternative 
called "time-reversal processing" has been proposed w^ith great success in acoustics. It has been shown that 
time-reversal is applicable to spatio-temporal phenomena that satisfy a wave-type equation and possess a time 
reversal invariance property. 

Dynamic focusing using time reversal is essentially a technique to "focus" on a reflective target or mass 
through a homogeneous or inhomogeneous medium that is excited by a broadband source. More formally, time- 
reversal focusing converts a divergent wave generated from a source into a convergent wave focused on that 
source (see Fig. 2). It can be thought of as an "optimal" spatio-temporal filter that adapts to «ie medium in which 
the wave front evolves and compensates for all geometric distortions white reducing the associated noise. The 
underlying theory and application of time-reversal techniques to acoustical problems have been developed along 
with a wide range of appiications and proof-in-principle experiments. These applications have yielded some 
exciting results in focusing through an inhomogeneous medium and offer an opportunity for many different 
applications. This approach has been demonstrated for the focusing and destruc^on of painful kidney stones in 
lithotripsy. Fortunately, unlike tissue mass, the stones are highly reflective and the most dominant scatterer in 
the kidney. 




Figure 2. Time Reversal Focusmg: After reception of scattered field, the temporal signals are reversed 
and retransmitted into the medium where the acoustic energy is focused on the mass. 

Model-Based Focusing 

An alternative to time reversal is the model-based approach that: 

• develops a model of the inhomogeneous medium including the mass under scrutiny from the results of 
quantitative imaging; 

• numerically propagates acoustic energy to the array from a virtual (fictitious) source located at the mass 
generating a set of synthesized multichannel time series; and 

• transmits the acoustic energy back into the medium to "focus" on the target mass. 

"Blind" time reversal that will focus on the strongest scattering mass in a completely unknown tissue medium 
without any a-priori information about the medium, mass or its location is clearly a risky endeavor. In contrast, 
the model-based approach uses the mode! of the medium (including the mass and its location) to synthesize the 
appropriate time series and focus at the correct location. Clearly, the major problem with this approach is the 



development of the appropriate model. We propose to apply quantitative imaging using tomographic 
reconstruction techniques to characterize the medium model arid art acoustic propagation algorithm to synthesize 
the required signals. We depict the model-based approach in Fig. 3. 



MODEL-BASED FOCUSING 




Figure 3. Model-Based Focusing: After quantitative imaging, the propagation model is characterized, 
temporal signals are generated, reversed and transmitted into the medium where the acoustic 
energy is focused on the mass. 



This completes the background of the invention of focusing techniques in tissue medium. 
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VIII. invention Description 

Description of the invention (you may also attach a paper). Piease include a sketch of the invention, if possible. 
Matched-Field Detection/Localization 



Matched-field processing (MFP) is considered by many to be an outgrowth of matched filtering in which a 
known signal such as a pulse in conventional ultrasound is transmitted into a medium and its return is to be 
detected from noisy measurements. Here a replicant of the pulse is convolved with the measurement to produce 
an optimal detection. When the pulse is unknown or cannot easily be measured or passive listening is assumed, 
then the replicant is no longer available and other methods must be used to generate the required replicant for 
c^ttmat detection. MFP uses a propagation model of Uie medium to generate the replicant for detection. It 
compares the model predicted field (replicant) propagated to the array position to the field actually measured at 
the sensor array to achieve the detection. In the localization problem, the MFP guesses at the position of a 
source, propagates it to the sensor array using the model and compares it to the measured field. That location 
with the maximum power is deemed the location of the source. A diagram of the MFP is shown in Fig. 3. After 
careful preprocessing to remove extraneous signals and noise, the data are ready for imaging, Each pixel in the 
Image representing a source or mass position is propagated to the sensor and its power or other feature is 
estimated to create the image. The threshold is applied to detect the presence of masses while their locations 
are determined by the coiresponding maxima. Thus, in this way MFP offers a reasonable approach to imaging 
for mass detection and localization, when a propagation model is available. MFP can be considered a 
generalization of standard (delay and sum) beam forming methods that replaces the plane or spherical wave 
models with more sophisticated propagation models. It is a general "model-based" methodology that can and 
has been extended to solve many different problems. 



MASS DETECnON/LOCAUZATION 




Fjgure 3. Mass Detection/LocalUation: Model-based approach using a 
matched-field processor. 



We begin our brief development of the processor with the overall field measured by a sensor or array of 
sensors arsd develop the basic signal models that will lead to a practical imaging technique. First, we develop 
the underlying mathematical relationships to characterize our measured wave field. 

Assume that the wave field resulting from the ultrasound satisfies the wave equation. The acoustic 
pressure at the -sensor is given by 

"{'y^')-0{r,,r^,t)*s{r,,t), (1) 

where 

ti{>y,t) is the ultrasonic wave field st the -sensor; G(r„/:j;?) is the Green's function of the medium at 
r(,r^ from the source-to-sensor at time t\ and .?():/,/) is the source at and time f. 

The actual sensor measurements are contaminated with gaussian random noise as well; therefore, we define 
the noisy sensor measurement field as 

z,{t)^u{r,j) + n,{t), {2} 
for the random noise contaminating the^ -th sensor. If we expand this expression over the entire L-etement 
sensor array, then we obtain the vgslsr measurement field 

z(t) = K(0+!i{f) =G(/)*s{r,,r) + «(/). (3) 
where HjGsC^"^ are the measurement, field signal, vrfiite gaussian noise vector of variance o-^i, Hie 
medium Green's function and the respective source (mass) terms. Using this generic measurement model 
representing the noisy wave field measured across the array, we next develop the matched-field (MF) processing 
approach. 

The underlying problem is to decide whether or not there exists a mass in the tissue specimen. Assume 
that we have the "known" replicant field signal, m(t) , generated from our developed model (discussed above). 
Our problem Is to detect a mass signal from the test specimen measurements. That is, we must solve the binary 
decision problem 

Ho:2(0 = s(') [noise only] 

H] : z(0 = OT(f)-*-«(0- [mass signal + noise] 
The solution to this problem is easily obtained from the Neyman-Pearson criterion and is given by the log- 
likelihood ratio test (LRT) 

A(^) = b Pr(zjHj)-ln Pr(zjHo) 5 In A , (5) 

Where Pr is the probability density function and X is the threshold of the test. This problem, assuming that the 
measurements are zero-mean, gaussian with variance o-fj leads to the decision function 

A{z) = -^^{,it)-m{t)) {z(t)-m{l))-z(t)z(t)j > In 1 . 

Expanding this expression and coHecting all data dependent terms, we obtain the sufficient statistic 

A [z) = m (02(0 i al In i+ 1 w' (0^(0 - • (6) 
Under the Neyman Pearson criterion, the threshold can be detemiined from the false alarm probability given by 

to a pre-selected value by solving for A, and X in Eq. 6. In the white, gaussian noise case, we have that 
Pr(A i //„)- n(o,c7^-i) which leads to the threshold lJoh93] 

A=yja-lEL<r^{Pp;,) (7) 
With the signal energy, E s ni\t)m{t) , cl> a unit variance gaussian distribution and L the number of sensors in 
the array. 

Note also that by a simple change of variables in f, it is easy to show that the sufficient statistic of Eq. 6 
is the well-known matched-filter solution with "matching" filter impulse response given in terms of our vector 
signal mode! of Eq. 6 by 



min ^>i(T -t), and A{z) ^ « 0-7 )*£(0 , (7) 
which is ^mply the time reversed, repficant of the known field. Recall also from matched-filter theory that the 
desired solution is to find the optimal filter at each sensor channel such that the output signal-to-noise ratio 
(SNR) is maximized, that is, the matched-filter is the solution (in time or frequency) to 



for (-^an appropriate inner product yielding again 



miO^uiT-t). (9) 
The important point here is that the matched-filter solution is simply the delayed, time reversed, replicant of 
the known field signal vector in the white, gaussian noise case. It is easy to extend this to the non-white noise 
case with the subsequent processor incorporating a pre-whitening filter (inverse of the noise covariance matrix) 
operation followed by the processor developed above. 

In our solution, we have assumed that the field vector, u({) , is completely known a priori. Suppose that the 
assumption is no longer true and we can characterize the unknown or missing parameters (e.g. amplitude, 
phase, etc.) by the embedded vector, 6, then our field vector becomes u(t;0) and therefore the "matching" 
vector is m{t;0}. The solution to this mass detection problem can be solved by composite hypothesis testing. 
In this case the test is 
IL : 2(0 = «ro 

H,:£(O = /«(/i0 + »(O ^ ' 

writh corresponding log-likelihood ratio 

A (2; £) = Iti Pr(£ I ^Jl, ) - In ?T(z I ^,Ho ) > In^. 

One solution to this problem is to estimate the parameter vector, 0 and then proceed as before which ieads to 
the generalized iog-likelihood ratio test (GLRT) 

max A {z; $) = max [in Pr(£ j ^,H, )] - max [in PrCg | |?,Ho )] > In . (11) 
Substituting m{i) in the previous relations, we have that 

A (2; 9) = m ((; e)zQ) ^' aj, in ^-^^m in 6)m{t;e) ^Xq. (12) 

The result implies that as we develop a solution to the mass detection problem, we must search over the 
unknown parameter set, {g} to maximize the log-likelihood uang the GLRT to "match" the model replicant field 
to the data measured across the sensor array. This approach then leads to matched-field detection. We search 
various parameter vectors and find that value d that leads to the maximum log-likelihood or equivalent 
maximum output SNR power defined by 

^^•^m-\r~-, ^- ^< ^- (13) 

Thus the detection of the mass is determined, when the set threshold is exceeded, if we assume (simply) that the 
mass can be represented by a spatio-temporal point source, then performing the prescribed convoiution with 
*0:.'f) = ^('-'i). we have that 

£(0 = G'(/)»<y(r-rJ^G'(/-rJ. (14) 



in terms of the matched-field approach, if we assume that the unknown parameters are the source or 
equivalently mass position, r^, then we see immediately that our matching or repiicant vector in the medium is 

given by &^ = r, - [x, x,] . position of the mass, thai is, the matched filter solution is 

!>Li'>i) = d{T-t + t,,i,). (15) 
Therefore, we can create output SNR "power" surface and detection scheme by forming the GLRT 

n«.x^(l) "i ^ (16) 

where F(i,)-^, \ ■ 

{miT-,d^)*m(T\e^)) (c? <^;^)*^(r;£,)\ 

Thus, the so-called "matched-field" detector/locaiizer uses an assumed position, 0, and the propagation model 

to produce the replicant, m{n0- The model replicant is then convoived (correlated) with the measurement, 

z(T) to produce the detection statistic, P{&^) which is compared to the threshold, to detect the presence of 

a mass at the pixel specified by the location parameter, 0. 

Iterative Time Reversal 

Iterative time-reversal techniques were discussed in a previous patent (see past references), here we briefly 
mention the methodology that will be used as an integral part of some of the subsequent inventions. In this 
section we discuss the time reversal approach to focusing. Recall thai time-reversal processing is a focusing 
technique that can be used to minimize the aberrations created by an inhomogeneous or random medium 
illuminated by propagating waves. This technique can be used to "focus" on the principal scatterer dominating a 
putse-echo response. The T/R technique simpfy processes the multichannel time series radiated from the region 
under Investigation, collects the array data, digitizes, time-reverses the temporal array signals and re-transmits 
them back through the medium to focus on each scatterer (see Fig. 2 above). 

In the decoupled scatterer case, i.e., each scatterer has a distinct (fixed) eigenvalue and eigenfunction 
associated with it, it is possible to perform the cycle "iteratively" by focusing on the strongest mass, receiving its 
scattered field and removing it from the time series data, then we could develop an iterative scheme as depicted 
in Fig. 3. Here the strongest scatterer signature is removed from the data sequentialiy until no other scatterers 
exist. One approach to this removal problem is suggested in the previous patent. 



ITERATIVE TIME REVERSAL FOCUSING 
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Figure 3. Iterative T/R focusing pitch-catch-remove'revetse sequence. 

Model-Based Focusing 

The model-based focusing approach; (1) develops a model of the inhomogeneous medium including the 
mass under scrutiny from the results of quantitative imaging; (2) backpropagates the iocaiized mass (source) to 
the array generating a set of synthesized array time series; and (3) transmits the time reversed acoustic energy 
back into the medium to "focus" on the target mass. In contrast to "blind" time reversal that will focus on the 
strongest scattering mass, the model-based approach uses the model of the medium (including the mass and its 
location) to synthesize the appropriate time series and focus at the correct location. We appiy quantitative 
imaging to characterize the medium model and an acoustic propagation algorithm to synthesize the required 
signals. We depicted the basic model-based focusing approach previously in Fig. 3. 

Interactive Model-Based Focusir)g 

Perhaps the simplest technique to localize a mass under scrutiny Is to enable the physician to examine 
the tissue image and select questionable regions for further more detailed investigations, just as a radiologist 
would do when examining x-rays for fractures. In this approach the physician uses, for example, an interactive 
light pen to select individual masses or zones requiring further detailed analysis. We depict this interactive 
model-based T/R focusing approach in the invention of Fig. 4, After selection of the mass, its position Is 
provided as input to the focus synthesizer that then generates the required time series from the fonward 
propagation/system mode!. After reversal the focusing signals are then transmitted into the medium and they 
coherently superpose at the de&red mass location for treatment. Conceptually, this approach is simple, but it 
relies heavily on the physician to select the appropriate masses for treatment or regions to be investigated more 
completely, 

Model-Based Iterative T/R Focusing 

The final approach is perhaps the most sophisticated in that it combines both the strength of the iterative 
T/R focusing and detection capability with the model-based focus synthesizer as shown in Fig. 5. Here we use 
the iterative time-reversal approach to "detect" the mass in a zonal region possibly selected by the physician 
(above). Once the mass is detected, it is iocaiized using the model-based, matched-field processor with the 
mode! developed from a quantitative image as before. After locaiization, the mass could be classified as benign 
or malignant (not part of this proposal), Once localized, the position of the mass is provided as input to the 
model-based focusing algorithm that produces the required set of time series. As before, the time series are 
reversed and transmitted into the medium to focus on the mass. After physical mass treatment, the procedure is 
repeated for the next mass to be treated. In principle, this approach appears to have the most appeal, since it 
employs the power of iterative time-reverser combined with the modei-t>ased focusing algorithms guaranteeing 
that the mass selected is to be treated. In any case, the algorithm discussed both model-based and time- 
reversal based offer the potential to perform noninvasive acoustic surgery. Before we close this section, it is 
Important to understand the basis of using the iterative T/R to detect the presence of the dominant mass, ttiat is, 
we must answer the question when has the algorithm converged to yield a detection. In the next section we 
address that new invention. 




Figure 4. Interactive model-based T/R focusing. A physician selects to region or zone to investigate and 
locates the mass under scrutiny providing mass position information to the focus synthesizer 
which generates the required time series that will be reversed and transmitted back into the 
tissue medium. 
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Figure S. Modet-Based T/R focusing: Iterative T/R detection and focusing is coupled with model-based iVIFP 
to detect, tocalize, and treat the mass under scrutiny by the physician. 



Iterative T/R Detection 

The development of a dominant mass detection algorithm using the T/R processor follows the same 
analysis as before using the iterative T/R models. In this section, we again develop a solution to the dominant 



mass (scatterer) detection prob. . Again we are assuming that the rece J field is contaminated by zero- 
mean, gaussiart noise of variance, (x^ , then the noisy array measurement becomes 

z(r:t)^R{r,t) + V{fj). (17) 
Our basic problem is to determine whether we have a single mass (scatterer) or equivalently has the 
iterative T/R processor "focused" on the dominant mass. If we assume this measurement model, then we must 
soive the following decision problem at each iteration, 
Hq ; Zf(r;0 = ^'/(r,-/) [Noise Only] 
H| : ^,{r,7) = R,{rQ;t) + Vi{r:t) [Signal + Noise] 
where z,,V„R, e R''*''''' with the array measurement for a single scatterer defined by 

Mh-;f)-^S,(r;i)*gM:'), (19) 
and 9,(ri,7) the scatterer return (scalar) associated with the /"^-iteration. Also, gi(r,7){S an iV^ -vector 
defined as the column of the A?,, xA'^ -Green's function matrix. This definition can be rewritten in expanded 
form as 



(18) 



fi(r;0-(?{r;/)*^(r;i)= g,(r,-0 g|(r;0 - g^^..,(r;f) 



(20) 



or performing these operations, we obtain 

«(r.-0 = rg„(r;0*^(«o:0+- + g;^,-i(r,-0*?(rw,-,;0]= t\kir;i)*q{r,U) (21) 

The solution to this problem is easily obtained from the Neyman-Pearson criterion as before in 5 given by the 
log-likelihood ratio test (LRT) 

A(z;) = ln Pr{z,(r;/)]H,)-ln Pr(^,(r;/)| Hq) > In Z , {22} 

Where Pr is the probability density function and X is the threshold of the test. This probiem, assuming that the 
measurements are contaminated by additive zero-mean, gaussian noise with variance a^l leads to the decision 
function 

A(2,) = ^^\^{z,(r;l)-R,{rj)) (z,(r;r)-/f,(r;0)-2;{r;0^/(r;f)] j In A . 

Expanding this expression and coiiecting all data dqjendent terms, we obtain the sufficient statistic 

(23) 

Under the Neyman Pearson criterion, the threshold can be determined from the false alarm probability. 

Note also that by a simple change of variables in f, it is easy to show that the sufficient statistic is the 
matched-filter solution with "matching" filter impulse response given in ternis of our vector signal model by 

/f;(r;r-/), and A{z,) = R,{r;(-T)*Ztij;i) , (24) 
which is simply the time reversed, replicant of the known field. The desired solution is to find the optimal filter at 
each sensor channel such that the output signal-to-noise ratio (SNR) is maximized, that Is, the matched-filter is 
the solution to 

- ~^-{Rlir,T)*R,(T-T)) ^{R;(r-OR,(T;4)d4) 

for Qan appropriate inner product. 



Applied to our problem, we see that the matching or replicant vector is given by, ;?,(«*o-'^-"0 .which is the time- 
reversed, received field induced by the dominant mass received at the array. Therefore, the detector of Eq. 25 
becomes 

1] ^ max SNR = '^ --i/— ^ A . (26) 

^ ^■■{R!ir,,T)*R,-{r„;r)) J 

The problem we have now is to estimate the required replicant. Ri%;(), in order to implement the 
optimal detector. We know that under certain conditions 

R^(r;t)^R,(Tg:l), for /-^iV,. 
where is the number of iterations required for the power method (T/R) to converge and is based on the ratio 
of the two largest scattering coefficients (eigenvafues). Thus, using the matched-filter theory [Joh93j developed 
above and the T/R focusing property, a pragmatic method of detection is to use the previous iterate, i?(-i(r;/) , 
produced during the "pitch-catch" sequence as the repli<^nt and continue the iteration until the output SNR does 
not change, that Is, 

^,O r;r-0.,(r.-0 \^ (27) 



i?,,2(r,T~0^M(r;f)J 



Clearly, P^^i as the T/R processor focuses on the strctngest mass, that is, |^-^jxiOO->iOO%. We 

demonstrate the performance of the detector on our homogenous medium simulation and show the sequence of 
convolutions during the convergence of the T/R to the dominant scatterer. Here we set the threshold, 
T = 99.5 % resulting in near perfect focusing and detection. Note that at each iteration the dominant mass 
return increases relative to the others. This completes the chapter on focusing, next we discuss the hardware 
and experiments available for this invention. 




Figure 6. T/R detection afgorithm performance on homogeneous medium with three masses: (a) Final 
(rigtit) and previous iterate, (b) Decision function (99.5 %}. (c) Received (raw) field at the array and T/R 
focused/detected mass. 

This approach of T/R detection can also be used to localize and map the mass as we discuss in our next part of 
the invention. 



Localization and Mass Mapping Invention 

jn this section we develop a tocaiization and mass detection technique (invention) based on the sdea of 
"wave front matching." Our approach is to first perform a homogeneous wave front match using a global 
technique to search for the best fit based on maximum power at a given location. The location (xy-position) 
output of this estimator then becomes the starting value for the local focusing algorithm that essentially performs 
a nonlinear least-squares fit over the region around the starting value. The focuser can be considered a zoom in 
approach to refine the grid and search. Our solution to this problem is shown in Fig. 7. Note that it is predicated 
on the fact that the T/R algorithm of the previous section has focused on the strongest scatterer and the 
decomposition algorithm has extracted it from the total received field data. Therefore our problem here is only to 
locate the position of this mass. 

Global Localization 

Our propagation model for this medium satisfies the homogeneous wave equation for a single scatterer, 
then under these assumptions the solution to the wave equation is that of a free space Green's function given by 



I r-r„ I , the Euclidean distance between the source at and the observation at r . 



(28) 








^ Focus j 













Fig ure 7. Mass Localization Algorithm using Global/Local iterations. 

Now returning to (28) using the homogeneous Green's function above and perfomriing the convolution, we obtain 
the wave field relation at the sensor as 



(29) 



where = 

If we now extend these models for a single scatterer at r„ obtained by the T/R processor over the A?; -element 
sensor array, we obtain the vector relations 

R{r,u) = g{r.lt)*s{r„U), (30) 



tf we choose to perform weighted delay-sum beam forming at the output of the array, then we obtain 



(31) 



Now if the beam former is steered to the correct scatterer location, then %=r„, W(,(0 = 4;irA/j, |r, "r<,|, and 
zg^to+ . The output is given by 

6/(r„;0-5(r^;0. . , 

and therefore, power output is maximized as 

P(r0) = \s{T^;i)f. (33) 
Thus, our approach to the gtofia/ search technique is based on matching the homogeneous wave front that is 
equivalent to performing delay-sum beam forming. Let us continue with our homogeneous example of the 
previous section and search over the dimensions of the part under evaluation by the following search technique: 

GLOBAL SEARCH ALGORITHM (HOMOGENEOUS WAVEFRONT) 



• decompose the part dimensions into pixels (ax^.a^^), i = j = l,—,N^ ; 

• for each (Ax„Ayy) calculate the corresponding time delay, r,(A) = !~^ ^""^ ' ■ Ajq, =/Ax. Avj =yAy. 

- I = ^/(^7^^Axf^(y, - JAyf ; 

• perform weighted sum-delay beam forming according to Eq. 4.6; 

• calculate the power, Pir^) , at the array output for each pixel; and 

• select the pixel of maximum power as the global search position estimate. 



We synthesized a point mass in a homogeneous medium of silica with sound speed 3.5 mm/usec under the 
same conditions of the previous example. We generated the field data as before with the true synthesized mass 
positioned at (12mm,6mm). The global search technique perfonns quite well (as expected) for the homogeneous 
case and the resulting power image is shown in Figure 15. Here we see the maximum located at approximately 
the true position. 

Local Focusing Approach 

Once we have a starting value resulting from the global search, we use Wiese estimates in a wave front 
matching algorithm. We set up the following nonlinear ieast-squares problem by first defining the error between 
the measured receiver array outputs, R(ru) , and the estimate, R{r;i) , that is, 

£(rgU) = RifU)-Rir,t) = R(r,t)- RitguJ), (34) 
which leads to the following cost function 

Using Eq. (28), we estimate the wave front received at the array by defining the following fonward 
propagation model, R(rJ). If we have a tomosieneoi/s model, then 
1 



(36) 



da(U) = \r-r,irJ)\ and r,(/J) = t-^^^^ for r,(ij) = ix,.,y,) . {Z7) 
The local focusing aSgorithm can be implemented by: 




X-posltton (mm) 

Figure 8. Global search approach to estimate mass (scatterei^ position using homogeneous wave front 
matching (true mass position (12mm,0mm}). 



LOCAL SEARCH ALGOmHM (HOMOGENEOUS CASE) 

• initialize the search with the initial global position estimates obtained from above, (i, J) (i,- , y, ) ; 

• estimate the con-esponding time delays, Tg{i,j) using (4.12) with ^rAr, =yAy, and 

• search over ail {ij}, i = l,...,Nj, J = using the polytope method [MAT93]; 

• esf/mafe for each {!,)> the mean-squared error (MSE), Jff{i,j} where sgii,j)sR(r,i)-Rij(T;tJ)] and 

• select the search portion estimate, j) == (jc* , y* ) correspondi ng to the minimum MSE. 



We used the same problem defined above and synthesized data at 3 dB SNR on a 32-elemetit array 
driven by a narrow pulse. The results of the combined global/locai localization iteration algorithm are show in 
Fig. 9 below. The estimate wave front at the estimated parameters is shown in 9a along with the mass map and 
the estimated position and mean-squared error function convergence {- 50 iterations). The results are quite 
good on this synthesized data set. This completes the localization algorithm based on the T/R processor. 

This completes the description of the overall invention. 
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Dr. Chambers is currently performing computer simulations of acoustic propagation in human 
tissue in support of a UCSB project. The continuation of this work is anticipated to taJke 25% 
time until the end of the project. 

Scope of Current LLNL Assignment 

This work will reduce time available for LLNL projects to 75%. However, Dr. Chambers 
responsibilities have included accepting work from a number of small projects and this new 
assignment would only limit his time to work on these. It would not affect his ability to work on 
ADAP T or other programs planned for FY03. 
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Date: Thu, 01 Aug 2002 15:33:49 -0700 

To: nelson20(ai lnl.gov 

From: Anna Wright <wright28®nnl .gov> 

Subject: Chambers, David - Charging 1990-09 

cc: dougan3@nnl.gov, chambers2@nnl.gov 

<x-fl owed>wayne: 

May David Chambers continue to charge against 1990-09 while we are catching 
up with the paperwork? UCSB (Megan Cron) has confirmed funds are available 
thru 11/15/02. I need to advise David. 

Thanks for the guidance. 



Anna Wright 
Employment Assistant 

Student and Faculty Employment Services 
Lawrence Livermore National Laboratory 
Phone: (925) 424-5480 
Fax: (925) 423-0894 
Emai 1 : wri ght28@l 1 n1 . gov 

</x-flowed> 
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to: Not chemistry chambers <chambers2@nn1 .gov> 

From: Anna Wright <wright28@nnl .gov> 

Subject: Chmbers, David - MLA Extension 

Cc : dougan3(&l 1 nl . gov 

<x-flowed>Dr. Chambers: 

We do need to complete paperwork to extend your Mul ti -Location appointment 
at UCSB. Per Megan cron, ucsb, the extension will be for the period 7/1/02 
- 11/15/02, when the funding expires. You will need to have your AD sign 
a memo extending your MLA appointment. Attached is a sample for you to 
use. YOU may use the period 8/1/02 - 11/15/02, since the original memo was 
for the period 8/1/02 - 7/31/02. I am coordinating the additional 
paperwork required CUPAY560 and amendment for the existing purchase order) 
with Megan Cron. 

Thank you for your cooperation. Please contact me or Molly Dougan, 3-8169, 

if you have questions. 

</x-flowed> 
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Memo-LLNL Home Ext. doc" 
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"James v. 'Jim' candy" <canclyl@llnl .gov> 
From: Pam Richmond <richmond2@llnl .gov> 
Subject: Fwd: MLA Questions 

<x-html> 

<html> 

<br> 

<blockquote type=cite class=cite cite>X~Sender : 
e228596®popgun . 1 1 nl . gov<br> 

X-Mailer: QUALCOMM Windows Eudora version 4.3.2<br> 

Date: Fri, 02 Aug 2002 15:52:16 -0700<br> 

To: Pam Richmond <richmond2@llnl .gov><br> 

From: "Amalia A. \"Molly\" Dougan" 

<dougan3(apopgun.llnl .gov&gt ;<br> 

Subject: mla Questions<br> 

Cc: "Anna M. wright" <wright28@llnl .gov><br> 
<br> 

Hi pam,<br> 
<br> 

Anna sent me your questions.   Here are my anwers below them.  

Let me know if you have any others or would like me to meet with you and 

the employees inquiring about MLA' s . --Mol ly<br> 

<br> 

<br> 

l.«  <b>Do the services performed for the other UC location have to 
be performed at the other uc location? could the services be performed at 
some other location, i.e., field location other than LLNL?<br> 
<br> 

</b>Yes, the services have to be performed at the UC location .  They 
cannot be performed at llnl or another field location unless it is one 
the UC location has set up for it's employees .  For example, uc 
Merced my set up a telecommuting site for employees to drive to and do 
work.  The campus has furnished the site, provided the equipment, 
etc.<br> 
<br> 

<b>2.  Who pays for travel and lodging expenses to the other UC 

location?<br> 

<br> 

</b>if the employee needs to be reimbursed for these expenses that's 
something they need to work out with the Host location.   UC travel 
policy does allow temporary employees to be reimbursed for these 
expenses. <br> 
<br> 

3.  <b>Can cost of travel and lodging be charged against the mla 

account?<br> 

<br> 

</b>No, only wages can be charged to the MLA account numbers. <br> 
<br> 

The important thing to remember is that an MLA is an assignment that is 
outside an employee's 'regular' job.  The Host Location is hiring 
our employee to do a job for them and paying for their expertise which 
may be similar to what they do here, but is not connected to the Lab in 
any way.  in other words it's not work for Others through a 
col 1 aborati ve agreement . </bl ockquote><br> 
</html> 
</x-html> 
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<x-flowed>Megan: 

Our Finance Department has requested authorization for Dr. David Chambers 
to continue charge time for MLA assignment services for UCS8. while we 
wait for the official purchase order amendment, can you email interim 
authorization for additional charges for the period 8/1/02 - 11/15/02 with 
the anotted doniar amount. There is currently a shortage in the amount 
to reimburse accumulated service charges, so Dr. Chambers is not able to 
continue to charge against this account without additional funding. 

Your assistance with expediting this matter is greatly appreciated. 

Thank you very much. 



Anna Wright 
Employment Assistant 

Student and Faculty Employment Services 
Lawrence Livermore National Laboratory 
Phone: (925) 424-5480 
Fax: (925) 423-0894 
Emai 1 : wri ght28@l 1 nl . gov 

</x-flowed> 
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To: nelson20(&ll nT.gov 

From: Anna Wright <wright28(ai Inl .gov> 

Subject: Chambers, David - UCSB Authorization 

Cc : dougan3(51 1 nl . gov , chambers2(ai 1 n1 . gov 

<x-fl owed>Wayne: 

Does this work for you? Please confirm. 



>X-sender : cron@mai 1 . ece . ucsb . edu 

>X-Mailer: QUALCOMM windows Eudora version 5.1 

>Date: Mon, 05 Aug 2002 13:24:18 -0700 

>To: Anna Wright <wright28@llnl .gov> 

>From: Meagan Cron <cron@ece.ucsb.edu> 

>Subject: Re: Chambers, David - mla Charge Authorization 

> 

>Hi Anna, 
> 

>This is to authorize spending for David Chambers salary for 
>7/l/02-10/31/02. There is currently $8950.00 left to be spent. This 
>should be appropriate since 25% of 8845.00 times four months is 
>$8845.00. I have excluded the month of November since there will not be 
>enough funds to cover him for that month. Please let me know if you need 
>any additional information. 
> 

>Thanks , 

>Meagan 

> 

> 

>At 11:00 AM 8/5/2002 -0700, you wrote: 

»Megan : 

» 

»Our Finance Department has requested authorization for Dr. David chambers 
»to continue charge time for MLA assignment services for UCSB. while we 
»wait for the official purchase order amendment, can you email interim 
»authorization for additional charges for the period 8/1/02 - 11/15/02 
»with the allotted dolllar amount. There is currently a shortage in the 
»amount to reimburse accumulated service charges, so Dr. chambers is 
»not able to continue to charge against this account without additional 
»funding. 
» 

»Your assistance with expediting this matter is greatly appreciated. 
» 

»Thank you very much. 

» 

» 

» 

»Anna Wright 
»Employment Assistant 

»Student and Faculty Employment Services 

»Lawrence Livermore National Laboratory 

»Phone: (925) 424-5480 

»Fax: C925) 423-0894 

»Emai 1 : wri ght28®l 1 nl . gov 

> 

>Meagan cron 

>Contracts & Grants/Financial Analyst 
>Electrical & Computer Engineering Department 
>university of California 
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chambers, David - UCSB Authorization.txt 
>santa Barbara, CA 93106 
>(805) 893-3939 



Anna Wright 
Employment Assistant 

Student and Faculty Employment Services 
Lawrence Livermore National Laboratory 
Phone: (925) 424-5480 
Fax: C925) 423-0894 
Emai 1 : wri ght28©l 1 nl . gov 

</x-flowed> 
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UCSB extension funding.txt 
X-Sender : e746690@popout . 11 n1 . gov 
X-Mailer: QUALCOMM windows Eudora version 5.0 
Date: Tue, 03 Sep 2002 10:41:25 -0700 
To: David Chambers <chamber52<§llnl .gov> 
From: Pam Richmond <richmond2®l1 nl .gov> 

Subject: Re: Fwd: Funding Agreement SB27222-0 was entered on 29-AUG-02 

<x-flowed>Thanks Dave. Actually our customer is UCSB, since the funding is coming 

to from UCSB not kci. That's why it's a mla and not a WFO. 
-Pam 

At 10:22 AM 9/3/02 -0700, you wrote: 

>Pam, 

> The dollar amount is consistent with what I expected to be sent 

> here, not UCSB. I don't know why they listed the customer as UCSB. It 

> should be KCI directly. I just returned from vacation and will call olsi 

> to confirm this, should they resend this with the customer changed? 
>-Dave 

> 

>At 10:32 AM 8/30/02 -0700, you wrote: 
»Dave, 

» I think this is the extension for your UCSB for KCI. Does the 

» dollar amount look right to you? 

»-pam 

» 

»>Date: Fri , 30 Aug 2002 09:45:28 -0700 (POT) 
>»From: Oracle Alert <appl mgrOai sq . 1 1 nl .gov> 
»>To: nelson20@llnl.gov, ri chmond2®llnl .gov 

»>Sub]ect: Funding Agreement SB27222-0 was entered on 29-AUG-02 
>» 

»>The funding agreement 5B27222-0 was entered in RWS on 29-AUG-02, and is now 

»>being processed for approval. 
»> 

»>Funding Agreement #: SB27222-0 
»>Agreement ID #: 1015731 
»>Proposa1 #: M9363 
»> 

»>customer: UC SANTA BARBARA 
»>5ponsor Amount: $27,450,00 
»> 

>»Expi ration Date: 
»> 

»>Resource Analyst: PAMELA RICHMOND 
»>Contract Administrator: WAYNE E NELSON 
»>Principal investigator: 

> ^ ^ 

>David H. Chambers 
>POB 808, L-154 
>Livermore, CA 945 51 
>925-423-8893 
>chambers2@nnl .gov 



</x-flowed> 
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Multi -Location Appointment.txt 
X-sender: e645974(apoptop.nnl .gov 
X-Mai1er: qualcomm windows Eudora version 5.0 
Date: Thu, 12 Sep 2002 11:17:27 -0700 
To: Anna Wright <wright28@llnl.gov> 
From: Wayne Nelson <nel 5on20(ai 1 nl . gov> 
Subject: Multi -Location Appointment 

Cc: richards3@llnl.gov, Not chemistry Chambers <chambers2®l Inl .gov> 

<x-htm1> 
<html> 

A purchase order has been received from HR for the following 

multi -location assi gnment : <br> 

<br> 

<b>Employee Name: <font co1or="#FF0000">Chambers , David<br> 
</font>Campus : UC/Santa Barbara<br> 
Account # Assigned: no change<br> 
Amount of P.O.: no change<br> 

Period Effective :<font color="#FFOOOO"> 7/1/02 through 10/31/02<br> 

</font>Resource Analyst: Dana Richard5<br> 

<br> 

</b>Please insure that all effort costs related to time worked at the 
campus for the period specified are either moved or charged to the 
account number assigned for this purpose .&nbsp ; The preferences and 
defaults can be adjusted in lite to include this account and the 
multi-loc earn type to make this easier on a weekly basis.   If you 
need help with this function in LITE please contact me and I would be 
happy to walk you through the process. Anbsp; <br> 
<br> 
<br> 

<x-si gsepxpx/x-si gsep> 

<font face= Coronet" si ze=6><b><i>Wayne Nelson <br> 

</i></font><font face="Tahoma" color="#808080">Lawrence Livermore 

National Laboratory <br> 

P.O. Box 808, L-435 <br> 

Livermore, CA  94550 <br> 

<br> 

</font><font face="Tahoma" size=2><i>Phone: 

</i></font><font face="Tahoma" size=2 color="#808080">C925)423-4443   
</font><font face="Tahoma" size=2><i>Fax: 

</i></font><font face="Tahoma" size==2 color="#808080">C925)422-3063</font><font 

face="Tahoma" size=2> 

</font></b><br> 

<font face="Coronet" size=6><b><i>Wayne Nelson <br> 

</i></font><font face="Tahoma" color="#808080">Lawrence Livermore 

National Laboratory <br> 

P.O. Box 808, L-435 <br> 

Livermore, CA  94550 <br> 

<br> 

</font><font face="Tahoma" size=2><i>Phone: 

</i></font><font face="Tahoma" size=2 color="#808080">(925)423-4443   
</font><font face="Tahoma" size=2><i>Fax: 

</i></font><font face="Tahoma" size=2 color="#808080">(925)422-3063</font><font 

face="Tahoma" size=2> 

</font></b></html> 

</x-html> 
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PATENT 

IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 



Applicant : 


James V. Candy et al. 


Docket No. : 


IL-10941 


Serial No. : 


10/661,249 


Art Unit : 


3737 


Filed : 


09/11/2003 


Examiner : 


James M. Kish 


For : 


DYNAMIC ACOUSTIC FOCUSING UTILIZING TIME REVERSAL 



DECLARATION UNDER 37 CFR §1.131 & §1.132 
Declaration by Eddie E. Scott of Prior Invention 
by TamevS V. Candy and David H. Chambers to Overcome Cite d I^eference 

Commissioner of Patents and Trademarks 
Alexandria, VA 22313-1450 

Dear Sir: 

I; hereby declare that: 

(1) . I am a citizen of the United States and a resident of Danville, 

California. 

(2) . I am employed by the Lawrence Livermore National Laboratory as a 
Assistant Laboratory Counsel and have been employed there from May 1, 1999 
to the present. 

(3) . I am empowered to act on behalf of the owner of the subject 

application, Lawrence Livermore National Laboratory, and I am the attorney 
representing the inventors James V. Candy and David H. Chambers (The 
Inventors), in the subject patent application. 

(4) . In the Office Action mailed October 2, 2007, the claims in the subject 
application were rejected over references that include the Kerbrat et al reference 
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(Transactions of Ultrasonics, August 2002). The date of the Kerbrat et al reference 
is on or about August 1, 2002. Accordingly, August 1, 2002 is the date that must 
be overcome. 

Overcoming the August 1, 2002 Date of Kerbrat et al Reference 

(5) , The Inventors made the invention described and claimed in the 
subject patent application (The Invention) in this country prior to the August 1, 
2002 date of the Kerbrat et al Reference. The Inventors made written descriptions 
of The Invention, The Inventors made tests of The Invention, and The Inventors 
disclosed Tlie Invention to others; all of the foregoing were done in this country 
prior to the August 1, 2002 date of the Kerbrat et al reference. Attached to the 
"Declaration by Inventors James V. Candy and David H. Chambers to Overcome 
Cited Reference" are documents supporting The Inventors statements. The 
documents support the fact that The Inventors made The Invention in this 
country prior to the August 1, 2002 date of the Kerbrat et al reference and that 
The Inventors were reasonably diligent in reducing their invention to practice 
during The Time Period. 

(6) The invention claimed in the subject patent application was described 
in a RECORD OF INVENTION document signed by James V. Candy and David 
H. Chambers, The Inventors of the subject application, on October 18, 2001. A 
copy of the RECORD OF INVENTION is attached. 

(7) . Upon information and belief. The Inventors and others were 
reasonably diligent in developing and reducing The Invention to practice in this 
country during the time period from August 1, 2002 date of the Kerbrat 
reference until September 12, 2002 when Applicants' Provisional Patent 
Application SN 60/410,575 was filed (hereinafter "The Time Period"). The 
""Declaration by Inventors James V. Candy and David H. Chambers to 
Overcome Cited Reference" and the attached documents support the statements 
that the Inventors and others were reasonably diligent in developing and 
patenting The Invention in this country during the time period from August 1, 
2002 when the Kerbrat reference was published until September 12, 2002 when 
The Inventors Provisional Patent Application was filed. 

Activity During the Time Period 

(8) . The Inventors and others continuously worked on testing, 
developing, and patenting The Invention during The Time Period from the 
Kerbrat et al reference date of August 1, 2002 until September 12, 2002 when 
Applicants' Provisional Patent Application SN 60/410,575 was filed. The attached 
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documents listed below support the statements that the Inventors and others 
continuously worked on testing, developing, and patenting The Invention 
during The Time Period. 

Screen Printout - Technology "Case No. IL10941" 
Screen Printout - "Technology /Event" Case No. IL 10941. 
September 6, 2002 emails regarding request for provisional 
September 12, 2002 Memo from Kathy Raymond to Nancy Stone 
September 12, 2002 Memo from Kathy Raymond to Rego & Rhodes 

(9) . The Inventors and others were reasonably diligent in developing and 
reducing The Invention to practice in this country during The Time Period from 
August 1, 2002 when the Kerbrat reference was published until September 12, 
2002 when Applicants' Provisional Patent Application SN 60/410,575 was filed. 

(10) Upon information and belief, the Industrial Partnership and 
Commercialization Office (IP AC) of the Lawrence Livermore National 
Laboratory continuously reviews inventions and prioritizes inventions for patent 
application filing. IPAC held monthly Invention Review Meetings during The 
Time Period and The Invention was reviewed at the Invention Review Meetings 
during The Time Period. The Invention was reviewed and prioritized by IPAC 
during The Time Period. The attached Screen Printout - Technology "Case No. 
IL10941" that IPAC held an invention review meeting during The Time Period 
and that the invention IL 10941 was reviewed at the meeting. The attached 
Screen Printout - "Technology /Event" Case No. IL 10941 shows that on August 
7, 2002 "R Kiefer - Paul Martin and I met with the inventor and we concluded 
that we would like to re-file the provisional. The inventors will be looking for 
funding. Campanies are still interested. There is an SBI proposal out." 

(11) . Upon information and belief, the Office of Laboratory Counsel 
(OLC) of the Lawrence Livermore National Laboratory held monthly Invention 
Review Meetings during The Time Period and The Invention was reviewed at 
the Meetings during The Time Period; 

(11). Upon information and belief, the Office of Laboratory Counsel 
(OLC) of the Lawrence Livermore National Laboratory held monthly meetings 
with the Industrial Partnership and Commercialization Office (IPAC) during The 
Time Period and The Invention was reviewed at the Meetings. 
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(13) . Upon information and belief, the Office of Laboratory Counsel 
(OLC) of the Lawrence Livermore National Laboratory, prepares patent 
applications for filing according to a priority list and the subject application was 
prepared by OLC covering The Invention during The Time Period according to 
the priority list. The attached documents "September 6, 2002 emails regarding 
request for provisional/' "September 12, 2002 Memo from Kathy Raymond to 
Nancy Stone/' and "September 12, 2002 Memo from Kathy Raymond to Rego & 
Rhodes" show that Applicants' Provisional Patent Application SN 60/410,575 
was being worked upon. 

(14) . The Inventors do not know and do not believe that the invention has 
been in public use or on sale in this country, or patented or described in a printed 
publication in this or any foreign country for more than one year prior to our 
application, and The Inventors have never abandoned their invention; 

(15) . I further declare that all statements made herein of my own 
knowledge are true and that all statements made on information and belief are 
believed to be true; and further that these statements were made with the 
knowledge that willful false statements and the like so made are punishable by 
fine or imprisonment, or both, under Section 1001 of Title 18 of the United States 
Code, and that such willful false statements may jeopardize the validity of the 
application or any patent issued thereon. 




January 35/ 2008 - (Signature) Declarant: Eddie E. Scott 
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RECORD OF INVENTION 



LLNL-I.RLG. 



This invention was made in the course of or under prime Contract No, \A/-7405-ENG"48 between the U.S. 
Department of Energy and the University of California, This Record of invention is prepared for the Office of the 
Assistant General Counsel for Patents, U.S. Department of Energy. 

I. Title of the Invention 

Dynamic Acoustic Focusing for Noninvasive Treatment 



LLNL File No, 



II. Inventor Information 



LLNL inventor(s) (F M L) 
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Hi. Abstract 

The ability to noninvasiveiy focus acoustical energy in tissue and directly on tissue masses (tumors, cysts, etc.) 
is the primary function of this invention. The objective is to provide the capability of focusing acoustic energy 
at a desired location for the purpose of treating tissue mass while minimizing the coilaterai damage in the 
surrounding tissue. This invention witt open new frontiers with the implication of noninvasive treatment of 
masses in the medical area along with the expanding technology of acoustic surgery. 



IV. Uses of the Invention 

List past uses, current uses and potential uses for your invention 

LLNL or Government uses or possibilities for use: 

Tissue mass removal, non-invasive tumor/cyst destruction, acoustic surgery, mass imaging, nondestructive 
evaluation of materials, secure communications, seismic detection of underground structures 



University of California 

Lawrence tivermore 
National Laboratory 
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V. Documents Describing the Invention 

Documents, publications, and presentations describing the Invention that you have pubiished or prepared for 
publication, or presented on the subject. Also inciude presentations and publications planned within one year 



Title/Subject 


Date 


Publication # 


Dynamic Acoustic Focusing for Noninvasive Treatment (Proprietary Doc.) 


8/9/01 


UCRL-PROP- 
145316 


Dynamic Focusing of Acoustic Energy for Nondestructive Evaluation 
(LORD Laboratory-Wide Final Report) 


11/01/01 


UCRL-JC- 
????? 















VI. Documents Describing Prior Art (Please include copies of these documents,) 

Related Documents (including patents, other pubiscations) Please include patent numbers, authors, title, 
publication date, etc. 



J. Candy, "Time Reversal Processing: An Approadi to the Scatterer Estimation Problem," UCRL-JC- 
124942,1996. 

J, Candy, "Dynamic Focusing of Acoustic Energy for Nondestructive Evaluation," LDRL-98-01, 1998. 

J. Candy, "Time Reversal Signal Processing: Background, Theory, and Applicaton," JASA-Voi. 101, No. 5, 
Pt. 2, p3089, 1997. 

Method and Apparatus for Dynamic Focusing of Ultrasound— IL-1 0,557 (pending No. WO 01/69283 A2) - 
James V. Candy, Sept, 20, 2001. 

J. Candy and D. Chambers, "The role of the time-reversal processor In acoustic signal processing," LLNL 
Report, UCRL-JC-141160. and J. Acoustical Soc. Amer., 108, 2483, 2000. 
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Vil. Background 

Background of the invention, including technical problems addressed by it: 

With the advent of high-speed digitizers, uttrafast computers, inexpensive memory, and the ability to 
construct dense acoustic arrays, the feasibility of noninvasive techniques of acoustic surgery offers a tantalizing 
alternative to current invasive techniques, in its simplest form, the focusing of acoustic energy to destructively 
treat a mass in surrounding tissue appears to be a reasonable approach to noninvasive surgery especiaHy if the 
medium is homogeneous and therefore can be characterized by attenuation and time delays. Thus, in this case, 
it is a simple matter to focus energy at a desired point in space. When the medium is inhomogeneous focusing 
at a desired focal point is more difficult unless some knowledge of the medium exists a-priori. There are two 
basic approaches that we discuss in this invention to focus energy in an inhomogeneous medium: (1) time 
reversal techniques, and (2) model-based techniques, and their combination. The time reversal approach 
requires no model in contrast to the model-based approach. We discuss both in this disclo^re. 

This invention discusses the capability of focusing acoustic energy at a desired location for ttie purpose 
of treating tissue mass while minimizing the collateral damage in the surrounding tissue. Our approach is 
summarized simply in f ig. 1. First, we must first detect the presence of a tissue mass applying acoustic energy 
propagated into the tissue using an array of ultrasonic transducers. The amount of energy scattered by the mass 
depends on its acoustic parameters (density, sound speed, attenuation, etc.). Once it is detected, the mass must 
be localized to determine Its position within the tissue medium. There are a variety of methods that can be used 
to perform this operation, but each make some underlying assumption about the characteristics of the medium 
(homogeneous, inhomogeneities with a homogeneous medium, etc.) leading to uncertainties, - Once detected 
and localized, temporal signatures must be developed to "drive" the array and focus increased energy back onto 
the mass through the medium. The increased energy generates heat, which essentially "cooks" the mass 
insuring its destruction. 

We demonstrate that confined focusing in tissue is pos^ble and discuss various approaches to achieve 
this focusing. Through "giobal" focusing, that is, inscMiifying a large region in the tissue the mass is detected and 
localized, "zonal" focusing is perfonned to extract or "zoom in" on the tissue mass under scrutiny. After it is 
decided to treat the mass, increased acoustic energy is transmitted back onto the mass. 




Figure 1. Conceptual Ultrasonic Focusing System for Noninvasive Mass Treatment: (1) Scatteiing 
tissue medium (3 large masses shown). (2) Focusing array. (3) Mass treatment via ultrasonic 
focusing. 



Time Reversal Focusing 



When a source propagates through a spatio-temporal medium, the resulting wave front is distorted, [f 
the medium is homoger^eous and the source resides in the near fieid, then a sphericai-type wave front evolves. 
But if the medium is inhomogeneous, then a distorted wave front results. In the first case, simple time-delay 
processing is sufficient to enhance the field at a given point; however, for inhomogeneous media the required 
time delays and amplitude are more difficult to estimate. The use of delay estimation and even adaptive delay 
estimation techniques become quite limited and unsuccessful in an inhomogeneous medium excited by a 
broadband incident field requiring an alternative s^proach to solve the focusing problem. A viable alternative 
called "time-reversal processing" has been proposed with great success in acoustics. It has been shown that 
time-reversal is applicable to spatio-temporal phenomena that satisfy a wave-type equation and possess a time 
reversal invartance property, 

Dynamic focusing using time reversal is essentially a technique to "focus" on a reflective target or mass 
through a homogeneous or inhomogeneous medium that is excited by a broadband source. More formally, time- 
reversal focusing converts a divergent wave generated from a source into a convergent wave focused on that 
source (see Fig. 2), It can be thought of as an "optimal" spatio-temporal filter that adapts to the medium in which 
the wave front evolves and compensates for ail geometric distortions white reducing the associated noise. The 
underlying theory and application of time-reversal techniques to acoustical problems have been developed along 
with a wide range of applications and proof-in-principle experiments. These applications have yielded some 
exciting results in focusing through an inhomogeneous medium and offer an opportunity for many different 
applications. This approach has been demonstrated for the focusing and destruction of painful kidney stones in 
lithotripsy. Fortunately, unlike tissue mass, the stones are highly reflective and the most dCMniinant scatterer in 
the kidney. 




Figure 2. Time Reversal Focusing: After reception of scattered field, the temporal signals are reversed 
and retransmitted into the medium where the acoustic energy is focused on the mass. 

Model-Based Focusing 

An alternative to time reversal is the model-based approach that: 

• develops a model of the inhomogeneous medium including the mass under scrutiny from the results of 
quantitative imaging; 

• numerically propagates acoustic energy to the array from a virtual (fictitious) source located at the mass 
generating a set of synthesized multichannel time series; and 

• transmits the acoustic energy back into the medium to "focus" on the target mass. 

"Blind" time reversal that will focus on the strongest scattering mass in a completely unknown tissue medium 
without any a-priorl information about the medium, mass or its location is clearly a risky endeavor. In contrast, 
the model-based approach uses the model of the medium (including the mass and its location) to synthesize the 
apprcH5riate time series and focus at the correct location. Clearly, the major problem with this approach is the 



development of the appropriate model. We propose to apply quantitative imaging using tomographic 
reconstruction techniques to characterize the medium model and an acoustic propagation algorithm to synthesize 
the required signals. We depict the modei-based approach in Fig. 3. 



MODEL-BASED FOCUSING 




Figure 3. Model-Based Focusing: After quantitative imaging, the propagation model is cliaracterized, 
temporal signals are generated, reversed and transmitted into the medium where the acoustic 
energy is focused on the mass. 

This completes the background of the invention of focusing techniques in tissue medium. 
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VIIL Invention Description 

Description of the invention {you may afso attach a paper). Please include a sl<etch of ttie invention, if possible. 



Matched-Field Detection/Localization 



lyatched-field processing (MFP) is considered by many to be an outgrowtii of matched filtering in which a 
known signal such as a pulse in conventional ultrasound Is transmitted into a medium and its return is to be 
detected from noisy measurements. Here a replicant of the pulse is convolved with the measurement to produce 
an optimal detection. When ttie pulse is unknown or cannot easily be measured or passive listening is assumed, 
then the replicant is no longer available and other methods must be used to generate the required replicant for 
optimal detection. MFP uses a propagation model of the medium to generate the replicant for detection, it 
compares the model predicted field (replicant) propagated to the array position to the field actually measured at 
the sensor array to achieve the detection. In the localization probiem, the MFP guesses at the position of a 
source, propagates it to the sensor array using the model and compares it to the measured field. That location 
with the maximum power is deemed the location of the source, A diagram of the MFP is shown in Fig. 3. After 
careful preprocessing to remove extraneous signals and noise, tiie data are ready for imaging. Each pixel in the 
image representing a source or mass position is propagated to the sensor and its power or other feature is 
estimated to create the image. The ttireshoid is applied to detect the presence of masses while their locations 
are determined by the coiresponding maxima. Thus, in this way MFP offers a reasonable approach to imaging 
for mass detection and localization, when a propagation model is available. MFP can be considered a 
generalization of standard (delay and sum) beam forming methods that replaces the plane or spherical wave 
models with more sophisticated propagation models, it is a general "modei-based" methodology that can and 
lias been extended to solve many different problems. 




Figure 3. Mass Det6ction/Loc0lization: Model-bssed approach using a 
matched-field processor. 



We begin our brief development of the processor with the overall field measured by a sensor or array of 
sensors ar^d develop the basic signal models that will lead to a practical imaging technique. First, we develop 
the underlying mathematical relationships to characterize our measured wave field. 

Assume that the wave field resulting from the ul^asound satisfies the wave equation. The acoustic 

pressure at the i^^ -sensor is given by 

u{ri,i) = G{r,.r,-,t)*s{r,,:), (1) 

where 

M (r^;/) is the ultrasonic wave field at the i"' -sensor; G (r.;,i:,;') is the Green's function of the medium at 
r(,r, from the source-to-sensor at time t; and j (r,;?) is the source at and time t. 

The actual sensor measurements are contaminated with gaussian random noise as weil; therefore, we define 
the noisy sensor measurement fieid a s 

z,{()^uir,u)^n,{,), (2) 
for the random noise contaminating the.f -th sensor. If we expand this expression over the entire iL-e!ement 
sensor array, then we obtain the vector measurement fieid 

£(0 - «(0 + "(0 = G(0 * ^ (r„ 0 + Hit). (3) 
where z,u, n.GeC^"' are the measurement, field signal, white gaussian noise vector of variance all, the 
medium Green's function and the respective source (mass) terms. Using this generic measurement mode! 
representing the noisy wave field measured across the array, we next develop the matched-field (Mf) processing 
approach. 

The underlying problem is to decide whether or not there exists a mass in the tissue specimen. Assume 
that we have the "known" replicant field signal, m(i) , generated from our developed model (discussed above). 
Our problem is to detect a mass signal from the test specimen measurements. That is, we must solve the binary 
decision problem 

Ha:^(0 = «(0 f noise only] 

H; : z(ii) = w(0 +«(')■ [mass signal + noisej 
The solution to this problem Is easily obtained from the Neyman-Pearson criterion and is given by the log- 
likelihood ratio test (LRT) 

A(z) = hi Pr{zlH5)-ln Pr(2lHo) > Inl, (5} 

where Pr is the probability density function and X is the threshold of the test. This problem, assuming that the 
measurements are zero-mean, gaussian with variance all ieads to the decision function 

Expanding this expression and collecting all data dependent temns, we obtain the sufficient statistic 

A (z) = m it)z(t) > cr^ In i (0«(f ) - A . (6) 
Under the Neyman Pearson criterion, the threshold can be determined from the false alarm probability given by 

to a pre-selected value by solving for A and A in Eq. 6. In the white, gaussian noise case, we have that 
Pr (-1 1 // J ~ N (o,a-„^l) which ieads to the threshold [Joh93] 

/1 = V^^'"'(^fa) (7) 
with the agnal energy, E s m^(t)rfi{t) , O a unit variance gaussian distribution and L the number of sensors in 
the airay. 

Note also that by a simple change of variables in t, it is easy to show that the sufficient statistic of Eq. 6 
is the well-known matched-fiiter solution with "matching" filter impulse response given in terms of our vector 
signal model of Eq. 6 by 



mii()^u(T-iX and A{i)=:u{i-T)*i{i), (7) 
Which is simply the time reversed, repficant of the known field. Recall also from matched-filter theory that the 
desired solution is to find the optimal filter at each sensor channel such that the output signal-to-noise ratio 
(SNR) is maximized, that is. the matched-filter is the solution (in time or frequency) to 

max — --=-4 ^ (8) 

for <|-)an appropriate inner product yielding again 

m{t)^u{T-t). (9) 

The important point here is that the matched-filter solution is simply the delayed, time reversed, replicant of 
the known field signal vector in the white, gaussian noise case. It is easy to extend this to the non-white noise 
case with the subsequent processor incorporating a pre-whitening filter (inverse of the noise covariance matrix) 
operation followed by the processor developed above. 

in our solution, we have assumed that the field vector, u{t) , is completely l<nown a priori. Suppose that the 
assumption is no longer true and we can characterize the unknown or missing parameters (e.g. ampiitude. 
phase, etc.) by the embedded vector, 6, then our field vector becomes u{V,&) and therefore the "matching" 
vector is n^V^d) . The solution to this mass detection problem can be solved by composite hypothesis testing. 
In this case the test is 

Hj :£(^) = ^nC/;^)^■^5^(0 ^ ' 

with corresponding log-likelihood ratio 

A(2;£) = [n Pr(£| ^,H,)-in ?T(z\eMo) ^ ln>^- 

One solution to this problem is to estimate the parameter vector, 0 and then proceed as before which leads to 
the generalized log-likeiihood ratio test (GLRT) 

maxA(£;^) = max[ln Pr(3 {£.H,)]-~max[k PrCgjgJlo)] >' Inl^. (11) 
Substituting jn{t;d)^m{0 the previous relations, we have that 

A(£;^)-m'C';0£(O i erf, ]n^-\-^i>l(ne)m{r,d)^h- (12) 
The result implies that as we develop a solution to the mass detection problem, we must search over the 
unknown parameter set, {$] to maximize the log-likelihood using the GLRT to "match" the model replicant field 
to the data measured across the sensor array. Tliis approach then leads to matched-field detection. We search 
various parameter vectors and find that value 0 that leads to the maximum log-likelihood or equivalent 
maximum output SNR power defined by 

maxP(f) = -V;-; ^ ^. (13) 

Thus the detection of the mass is determined, when the set threshold is exceeded. If we assume (simply) that the 
mass can be represented by a spatio-temporal point source, then performing the prescribed convolution with 
J (i'.'f )'=^ii-'s)> we have that 

zii)^G{i)*S{t-l,)^G{l~t,). (14) 



In terms of the matched-field approach, if we assume that the unknown parameters are the source or 
equtvalently mass position, r,, then we see immediately that our matching or replicant vector in the medium is 

given by -r^ ~[x, y^] , the position of the mass, that is, the matched filter solution is 

^:0;^)-G(r-r (15) 
Therefore, we can create output SNR "power" surface and detection scheme by forming the GLRT 

n...P{is) "i 4 (16) 



where 



Thus, the so-cailed "matched-field" detector/localizer uses an assumed position, 0, and the propagation model 
to produce the replicant, ni{t,&). The model replicant is then convolved (correiated) with the measurement, 
z{T) to produce the detection statistic, P{e^) which is compared to the threshold, to detect the presence of 
a mass at the pixel specified by the location parameter, d . 

Iterative Time Reversal 

Iterative time-reversal techniques were discussed in a previous patent (see past references), here we briefly 
mention the methodology that will be used as an integral part of some of the subsequent inventions. In this 
section we discuss the time reversal approach to focusing. Recall that time-reversal processing is a focusing 
technique that can be used to minimize the aberrations created by an inhomogeneous or random medium 
illuminated by propagating waves. This technique can be used to "focus" on the principal scatterer dominating a 
pulse-echo response. The T/R technique simpiy processes the multichannel time series radiated from the region 
under investigation, collects the array data, digitizes, time-reverses the temporal array signals and re-transmlts 
them back through the medium to focus on each scatterer (see Fig. 2 above). 

in the decoupled scatterer case, i.e., each scatterer has a distinct (fixed) eigenvalue and eigenfunction 
associated with it, it is possible to perform the cycle "iterativeiy" by focusing on the strongest mass, receiving its 
scattered field and removing it from the time series data, then we could develop an iterative scheme as depicted 
in Fig. 3. Here the strongest scatterer signature is removed from the data sequentially unti! no other scatterers 
exist. One approach to this removal problem is suggested in the previous patent. 



ITERATIVE TIME REVERSAL FOCUSING 




Figure 3. Iterative T/R focusing pitch'catch'remove-reverse sequence. 

Model-Based Focusing 

The model-based focusing approach; (1) develops a model of the inhomogeneous medium including the 
mass under scrutiny from the results of quar^titative imaging; (2) backpropagates the locaiized mass (source) to 
the array generating a set of synthesized array time series; and (3) transmits the time reversed acoustic energy 
back into the medium to "focus" on the target mass, in contrast to "blind" time reversal that will focus on tlie 
strongest scattering mass, the model-based approach uses the model of the medium (including the mass and its 
location) to synthesize the appropriate time series and focus at the correct location. We apply quantitative 
imaging to characterize the medium model and an acoustic propagation algorithm to synthesize the required 
signals. We depicted the basic model-based focusing approach previously in Fig. 3, 

Interactive Model-Based Focusing 

Perhaps the simplest technique to localize a mass under scrutiny is to enable the physician to examine 
the tissue image and select questionable regions for further more detailed investigations, just as a radiologist 
would do when examining x-rays for fractures, in this approach the physician uses, for example, an interactive 
light pen to select individual masses or zones requiring further detailed analysis. We depict this interactive 
model-based T/R focusing approach in the invention of Fig. 4. After selection of the mass, its position is 
provided as input to the focus synthesizer that then generates the required time series from the forward 
propagation/system model. After reversal the focusing signals are then transmitted into the medium and they 
coherently superpose at the desired mass location for treatment. Conceptually, this approach is simple, but it 
relies Heavily on the physician to select the appropriate masses for treatment or regions to be investigated more 
completely. 

Model-Based Iterative T/R Focusing 

The final approach is perhaps the most sophisticated in that it combines both the strength of the iterative 
T/R focusing and detection capability with the modei-based focus synthesizer as shown in Fig. 5. Here we use 
the iterative time-reversal approach to "detect" the mass in a zonal region possibly selected by the physician 
(above). Once the mass is detected, it is localized using the modei-based, matched-fieid processor with the 
model developed from a quantitative image as before, After localization, the mass could be classified as benign 
or malignant (not part of this proposal). Once localized, the position of the mass is provided as input to the 
model-based focusing algorithm that produces the required set of time series. As before, the time series are 
reversed and transmitted into the medium to focus on the mass. After physical mass treatment, the procedure is 
repeated for the next mass to be treated. In principle, this approach appears to have the most appeal, since it 
employs the power of iterative time-reverser combined with the modei-based focusing algorithms guaranteeing 
. that the mass selected is to be treated. In any case, the algorithm discussed both modei-based and time- 
reversal based offer the potential to perform noninvasive acoustic surgery. Before we close this section, it is 
important to understand the basis of using the iterative T/R to detect the presence of the dominant mass, that is, 
we must answer the question when has the algorithm converged to yield a detection. In the next section we 
address that new invention. 



INTERACTIVE MODEL-BASED T/R FOCUSING 




Figure 4. Interactive model-based T/R focusing. A physician selects to region or zone to investigate and 
locates the mass under scrutiny providing mass position information to the focus synthesizer 
which generates the required time series that wili be reversed and transmitted bacl< into the 
tissue medium. 
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Figure S. Model-Based T/R focusing: Iterative T/R detection and focusing is coupled with model-based MFP 
to detect, localize, and treat the nnass under scrutiny by the physician. 



Iterative T/R Detection 

The development of a dominant mass detection algorithm using the T/R processor follows the same 
analysis as before using the iterative T/R models, !n this section, we again develop a solution to the dominant 



mass (scatterer) detection prob. . Again we are assuming that the rece J fieid is contaminated by zero- 
mean, gaussian noise of variance, , then the noisy array measurement becomes 

z(r,0-/?(r,r) + K(r;0. (17) 
Our basic problem is to determine whether we have a single mass (scatterer) or equivalently has the 
iterative T/R processor "focused" on the dominant mass. If we assume this measurement model, then we must 
sofve the following decision problem at each iteration, 
Ho : z;(r;0 = y,{r;0 [Noise Only] 
H, : z,ir;t) = «,(ro;/) + )^(r,-;) [Signal + Noise] 
where Zi,Vj,R, e R''^'''' with the array measurement for a single scatterer defined t>y 

&(rA..-/)^gt(r;0*9;{r*.-0, (19) 
and (j^Cr^.r) the /c"' scatterer return (scalar) associated with the /"-iteration. Also, gi(r,7)is an A^^ -vector 
defined as the /c"" column of the //f xJV, -Green's function matrix. This definition can be rewritten in expanded 
form as 



(18) 



/?(r;/) = G(r,0*9(r;0 = [g.(r;0 g,(r,-0 - g,- ..,(r;f) J 



?(ro,-?) 



(20) 



or performing these operations, we obtain 

/?(r;() = [g,(r;/)*^(ro;0 + - + g;/,-i(«'.0*?(r;/,-i;0]= f gi(i*,O*?(r*;0 (21) 

The solution to this problem is easily obtained from the Neyman-Pearson criterion as before in 6 given by the 
log-likelihood ratio test (LRT) 

A(2;)-ln Pr(2,(r;01H,)-!n Pr(z,(r;0|Ho) > in X , (22) 

Where Pr is the probability density function and X is the threshold of the test. This problem, assuming that the 
measurements are contaminated by additive zero-mean, gaussian noise with variance a^l leads to the decision 
function 

= --l^[(z,(r;/)-«,(r;/))' {z,(Tj)-Ri(rj))- z;(r;/)7,(r;;)] i In i - 
Expanding this expression and collecting all data dependent terms, we obtain the sufficient statistic 

A(zi) = zUr'ORiir-i) "> a^. In i+i/?;(r;/)i?,(r:.0 = X , (23) 

Under the Keyman Pearson alterion, the threshold can be determined from the false alarm probability. 

Note also that by a simple change of variables in f, it is easy to show that Uie sufficient statistic is the 
matched-filter solution with "matching" filter impulse response given in terms of our vector signal model by 

Ri{r\T-l),mA A{z,)^ R,{r;t -T)*z,{r;l) , (24) 
which is simply the time reversed, replicant of the known fieid. The desired solution is to find the optimal filter at 
each sensor channel such that the output signal-to-noise ratio (SNR) is maximized, that is, the matched-filter is 
the solution to 

.ax.^.. f(-^n-.(r.r)r ^^^^ 



for (•)an appropriate inner product. 



Applied to our problem, we see that the matching or replicant vector is given by, i?,(r(,;7--0,which is the time- 
reversed, received fieid induced by the dominant mass received at the array. Therefore, the detector of Eq. 25 
becomes 



R =s max SNR - 



«(r,;7-)*z,(r;r)) > ^ ^^6] 



^-{R;ir,;T)*Ri(r„;T)) 

The problem we have n<m is to estimate the required replicant, «i(ro;f). fwder implement the 
optimal detector. We know that under certain conditions 

R,{r:t)=>R,ir^;ll for i-^N,, 
where A^, is the number of iterations required for the power method {T/R) to converge and is based on the ratio 
of the two largest scattering coefficients (eigenvalues). Thus, using the matched-filter theory [Joh93! developed 
above and the T/R focusing property, a pragmatic method of detection is to use the previous iterate, i?i_|(r,-/) , 
produced during the "pitch-catch" sequence as the replicant and continue the iteration until the output SNR does 
not change, that is, 

( pA ( i;. „(r;y-/K(r,-0 ^,.^ (27) 

Clearly, /^_! as the T/R processor focuses on the strongest mass, that is, We 

demonstrate the performance of the detector on our homogenous medium simulation and show the sequence of 
convolutions during the convergence of the T/R to the dominant scatterer. Here we set the threshold. 
j = 99_5 0/^ resulting in near perfect focusing and detection. Note that at each iteration the dominant mass 
return increases reiative to the others. This completes the chapter on focusing, next we discuss the hardware 
and experiments available for this invention. 



Figure 6. T/R detection algorithm performance on homogeneous medium with three masses: (a) Final 
(right) and previous iterate, (b) Decision function (99,5 %}. (c) Received (raw) field at the array and T/R 
focused/detected mass. 

This approach of T/R detection can also be used to localize and map the mass as we discuss in our next part of 
the invention. 



Locaiization and Mass Mapping Invention , u ^ .t,^ r 

m this section we develop a tocaiizalion and mass detection technique (invention) based on the idea ot 
"wave front matching." Our approach is to first perform a homogeneous wave front match using a global 
technique to search for the best fit based on maximum power at a given location. The location (xy-position) 
output of this estimator then becomes the starting value for the locel focusing algorithm that essentially performs 
a nonlinear least-squares fit over the region around the starting value. The focuser can be considered a zoom m 
approach to refine the grid and search. Our solution to this problem is shown in Fig. 7. Note that it is predicated 
on the fact that the T/R algorithm of the previous section has focused on the strongest scatterer and the 
decomposition algorithm has extracted it from the total received field data. Therefore our problem here ts only to 
locate the position of this mass. 

Global Localization , , . , ^ 

Our propagation model for this medium satisfies the homogeneous wave equation for a single scatterer, 
then under these assumptions the solution to the wave equation is that of a free space Green's function given by 



(28) 



with I r~ rj , the Euclidean distance between the source at and the observation at r . 
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Figure 7. Mass Localization Algorithm using Global/Local iterations. 

Now returning to (28) using the homogeneous Green's function above and performing the convolution, we obtain 
the wave field relation at the sensor as 



where r,=l^i::>i. 

V 

If we now extend these models for a single scatterer at r,, obtained by the T/R processor over the /V^ -element 
sensor array, we obtain the vector relations 

J?(r„;/) = g{r,;/)*5(r,;0, (30) 



where gifo'>') = 



4jr|r, -r^i 



If we choose to perform weighted deiay-sum beam forming at the output of the array, then we obtain 
bm^O^-^rtMO m>'-'o-K,+ r^)- (31) 

Now if the beam former is steered to the correct scatterer location, then Vg=r^. w5(£) = 4;rA/i,|r^ -r^| , and 

tg=to+ T, . The output is given by 

^,/{r,;O = ^(r.;0. . ^^^^ 

and therefore, power output is maximized as 

Thus our approach to the global search technique is based on matching the hwnogeneous wave front that is 
equivalent to performing delay-sum beam forming. Let us continue with our homogeneous example of the 
previous section and search over the dimensions of the part under evaluation by the following search technique: 

GLOBAL SEARCH ALGORITHM (HOMOGENBOUS WAVEFRONT) 



• decompose the part dimensions into pixels [isx„Ayj), i = J = i,-, A 

• for each {&x,Ayj) calculate the corresponding time delay, 7r,{A) = ^'^ 

• perfomi weighted sum-delay beam forming according to Eq. 4.6; 

• calculate the power, P{r^ ) , at the array output for each pixei; and 

• select the pixel of maximum power as the global search position estimate. 



Ax;- = iAs, Ay J = J&y, and 



We synthesized a point mass in a homogeneous medium of silica with sound speed 3.5 mm/usec under the 
same conditions of the previous example. We generated the field data as before with the true synthesized mass 
positioned at (12mm, 6mm), The global search technique performs quite well (as expected) for the homogeneous 
case and the resulting power image is shown in Figure 15. Here we see the maximum located at approximately 
the true position. 

Local Focusing Approach , ■ ^ * 

Once we have a starting value resulting from the giobai search, we use these estimates )n a wave front 

matching algorithm. We set up the following nonlinear least-squares problem by first defining the error between 

the measured receiver array outputs, R(r,t) , and the estimate, Rir;t) , that is, 

£(ry,t) » Rir;t)-kr\t) = /?(r;/)- ^(r;,;/,^) . (34) 

which leads to the following cost function 

Using Eq. (28), we estimate the wave front received at the array by defining the following fonward 
propagation model, i?(r;/) . If we have a homogeneous model, then 



Rirj,0)= j?(r;/-T,(/J)), 



(36) 



do{iJ)--=\r-r,{i,j)\ and T,{iJ)J^^^^^ for r,(/,7) = • (37) 

The local focusing algorithm can be implemented by: 
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Figure 8. Global search approach to estimate mass (scatterer) position using homogeneous wave front 
matching (true mass position {12mm,6mm}). 



LOCAL SEARCH ALGORfTHM (HOMOGENEOUS CASE) 

• initiaiize the searcli with the initial global position estimates obtained from above. r^(i,y) ^ (x,, j;,) ; 

• estimate the corresponding time delays, vgiij) using (4.12) m\h =i&x, yj ='J&y, and 

• search over ati {ij}, ( = !,...,a'^,7:=i,...,#^ using the poiytope method [MAT93]; 

• esf//nafe for each the mean-squared error (WISE), JgiiJ) where sg(i,j)^Riryi)~Rij(r;i,d); and 

• se/ecf the search position estimate, fo(/,j) = {x*,>'*) corresponding to the minimum MSE. 

We used the same problem defined above and synthesized data at 3 dB SNR on a 32-element array 
driven by a narrow pulse. The results of the combined gtobal/tocal localization iteration algorithm are shown in 
Fig, 8 below. The estimate wave front at the estimated parameters is shown in 9a along witii the mass map and 
the estimated position and mean-squared error function convergence (~ 50 iterations). The results are quite 
good on this synthesized data set. This completes the localization algorithm based on the T/R processor. 

This completes the description of the overall invention. 
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Figure 9. Gtobal/Locat Localization Algorithm Wave front Matching at MB SNR. 
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Date: Fri, 06 Sep 2002 10:53:37 -0700 

To: thompson52@i!nl.gov, stone26@tlnl.gov, kraymond@llnl.gov 

From: Bert Weis <weis1 @llnl.gov> 

Subject: IL-10941, Dynamic Acoustic Focusing for Non-invasive 
Treatment, Candy et al, Request for a Provtsionai. 
Cc: rhodes3@llnl.gov, berson1@llnl.gov, martin22@llnl,gov, 
Connie Pitcock <pitcock1 @llnl.gov> 

HiAl: 

This is to request the preparation and filing of a Provisional Patent Application on the above case. 

There was an earlier Provisional which was filed on 1 1/8/01 . it has been decided not to perfect 
this Provisional principaliy because no funding could be secured for the reduction to practice and 
further development and/or validation of the invention. However, the concepts of the invention are 
still considered to be viable and it may stili be possible to obtain funding. It is therefore desired to 
preserve as early a priority date as possible and on the other hand delay filing unti! such further 
work has been done. We therefor request the filing of another Provisional. 

The invention has not yet been published 

Thanks 
Bert 



Printed for Kathy Raymond <raymond3@llnl.gov> 



1 



University of California 
Lawrence Livermore National Laboratory 
Intellectual Property Law Group 

September 12, 2002 

Mail Station: L-703 

Extension: 3-9034 



Interoffice Memorandum 

TO: Janet Rego L-376 

Nina Rhodes L-795 

FROM: Kathy Raymond 

SUBJECT: Provisional Application IL~10941 (Second Filing) 



Enclosed for your records are copies of the Provisional Patent Application, the 
Provisional Application for Patent Cover Sheet, and our return postcard. 

Nina; Pursuant to Jan Tulk's policy dated January 7, 1998, concerning 
provisional and bar date cases, this memo serves as notice of our filing a 
provisional. 



Enclosures: 
as noted 



University of California 
Lawrence Livermore National Laboratory 
Intellectual Property Law Group 

September 12, 2002 

Mail Station: L-703 

Extension: 3-9034 

Interoffice Memorandum 
TO: Nancy Stone 

FROM: Kathy Raymond 

SUBJECT: USE OF PATENT GROUP DEPOSIT ACCOUNT (12-0695) 

As of this date/ 1 have charged the following amount(s) to the Patent Group deposit 
account: 

Amount IL lS f umber Attorney Type of Action 

$80.00 IL-10941 EES Provisional Application 



